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For dissipation-free photon-photon interaction at the single photon level, we analyze one-photon transi- 
tion and two-photon transition induced by photon pairs in three-level atoms using two-photon wavefunc- 
tions. We show that the two-photon absorption can be substantially enhanced by adjusting the time corre- 
lation of photon pairs. We study two typical cases: Gaussian wavefunction and rectangular wavefunction. 
In the latter, we find that under special conditions one-photon transition is completely suppressed while 
the high probability of two-photon transition is maintained. 

KEYWORDS: two-photon absorption, photon pair, quantum optics, two-photon wavefunction, time cor- 
relation, entanglement-induced two-photon transparency 
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1. Introduction 

Two-photon absorption is one of the most fundamental 
nonlinear processes, and it sometimes reveals the quantum 
nature of light. 1 If a two-photon transition is achieved at a 
single-photon level, it can be possible to control the fate of a 
single photon by controlling the presence or absence of the 
other single photon. Such photon-photon interaction may be 
applicable to quantum information technologies, such as the 
development of optical switches for two photons and Bell- 
state analyzers 2 among others. 3-8 In these applications, the 
medium must be such that it absorbs the two photons and 
not one of them. In three-level systems having a ground 
state, intermediate state, and an excited state, two-photon ab- 
sorption and one-photon absorption to the intermediate level 
can occur simultaneously. The two-photon absorption can 
be enhanced by decreasing the detuning to the intermediate 
state; however, this may also result in an increase in the one- 
photon transitions to the intermediate level. Several methods 
for enhancing the two-photon absorption while suppressing 
the one-photon absorption have been proposed. The use of 
electromagnetically induced transparency (EIT) is one ap- 
proach for realizing the suppression of one-photon transition 
while maintaining the strong nonlinearity in four-level atoms 
using auxiliary light called the coupling light. 911 The cav- 
ity QED is another promising method for obtaining effective 
two-photon absorption. 12 

It is also possible to enhance two-photon absorption by 
tailoring the quantum state of photons; 13, 14 this does not in- 
volve the use of any external apparatus such as a cavity or 
a control light, as in the case of the EIT system. To ob- 
tain efficient two-photon excitations, the two photons must 
satisfy the following two conditions: they should be close 
to each other, and the total linewidth of the two photons 
(two-photon linewidth) must be less than the linewidth of 
the excited state. Satisfying the first condition enables the 
occurrence of instantaneous transition via the virtual state, 
and the second one ensures that two-photon resonance oc- 
curs for a long time. An ultrashort light pulse does not sat- 
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isfy the second condition because of its broad spectrum. In 
contrast, twin photons produced by spontaneous paramet- 
ric down-conversion (SPDC) processes using a continuous 
pump laser satisfy both the conditions and can effectively 
undergo two-photon transitions. However, the suppression 
of the one-photon transitions to the intermediate state and 
enhancement of the two-photon transitions have not been in- 
vestigated in detail. 

In order to evaluate the transition probabilities of photon 
pairs, full quantum treatment of the system is required. We 
consider a propagating light beam consisting of continuous 
modes, which requires a multimode description. To carry 
out the full quantum multimode analysis of this light beam, 
we introduce two-photon wavefunctions that can represent 
arbitrary two-photon states. 15-19 Our formulation allows us 
to analytically estimate not only two-photon absorption but 
also the one-photon absorption. 

In Sec. 2, we describe the formulation both the two-photon 
and one-photon absorption for arbitrary two-photon states 
using two-photon wavefunction. In Sec. 3, we compare the 
two-photon and one-photon absorptions for two types of 
time-correlated photon pairs. Finally, we show that by ad- 
justing the correlation time, the single-photon loss can be 
suppressed without decreasing the high probability of two- 
photon absorption. 

2. Theory 

2.1 Field operators in free space 

In order to investigate the interaction of atoms with an 
electromagnetic field propagating in one dimension, we de- 
fine the electric field in continuous modes as 



E(z,t) 



ha; 

do;,, 



a(u>) e 



-iaj(i — z J c) 



H.c. 



(1) 



where A denotes the cross section of the beam. 15 The anni- 
hilation operator a(u>) satisfies the following commutation 
relation: 



[a(w), a* («</)] =8{oj-u'). 
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Fig. 1. Level diagrams for (a) one-photon transition and (b) two-photon 
transition. 



If the bandwidth of the field is assumed to be considerably 
less than the carrier frequency Hi, we can simplify eq. (1) as 



E(z,t) = £a(t- z/c) +H.c, 



(3) 



where 
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The Fourier transformed operator a(t) has the same commu 
tation relation as a(iv): 

[a{t),tf{t')]=6{t-t'). 



(6) 



2.2 Single-photon transition 

Before discussing the two-photon transitions, we first con- 
sider the interaction between a wavepacket containing a sin- 
gle photon and a two-level atom, as shown in Fig. 1(a). The 
initial state of the single-photon wavepacket is represented 
as and the atom at z = is prepared in the ground state 
\g). We introduce a simplified notation \g,ip) = \g)\ip) to 
represent the composite system. The Hamiltonian in the in- 
teraction picture is 



H(t)=p{t)E(t), 



where 



p(t)=V\e}(g\e int +H.c. 



(7) 



(8) 



represents an atomic dipole oscillating at the transition fre- 
quency fi and the electric-dipole transition matrix element 
given as V = —e(e\r\g). 16 

By applying the first-order perturbation theory, the proba- 
bility amplitude for the atom to exist in the excited state \e) 
after the passage of the wavepacket is given as 



ai 



1 r 



dt(e,0\H(t)\g,iP), 



(9) 



where |e, 0) = |e)|0) with the vacuum state |0). We assume 
that the transit time of the wavepacket is shorter than the re- 
laxation time of the excited state |e). On substituting eqs. (3) 
and (7) in eq. (9), we obtain 
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In this derivation, we use the rotating wave approximation. 
We introduce the one-photon amplitude 



i;{t) = (0|a(i)|V), 



(11) 



whose square is proportional to the probability of photon de- 
tection at time t. 16 Its Fourier transform is expressed as 

i r°° 

*( W ) = — / dt^tje"** (12) 

the square of which represents the spectral intensity of the 
wavepacket; then, the probability of excitation, Pi = |ai| 2 , 
can be expressed as 

P 1 =2nr\^(n)\ 2 1 (13) 

where r = £ 2 V 2 /h 2 . Thus, the probability of the one- 
photon absorption is determined by the spectral component 
| 2 at the atomic transition frequency of ft. 

2.3 Two-photon transition 

In this section, we discuss the two-photon transitions in- 
duced by a pair of photons, as shown in Fig. 1 (b). Let us 
consider that photon 1 in a mode induces the lower transi- 
tion, \g) — ► | a) while photon 2 in another mode is responsi- 
ble for the upper transition, |a) — > \b). 

The interaction Hamiltonian for this scenario is repre- 
sented as 



where 



H{t)=p 1 (t)E 1 (t)+p 2 (t)E 2 (t), 
p 1 (t)=V a \a)(g\e in « t 
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To estimate the two-photon excitation probability, we use the 
second-order perturbation theory. The second-order compo- 
nent of the time evolution operator is expressed as 
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dt 2 dtifl-(t 2 )#(ti)0(t2-ti), (17) 



where 9(t) is the Heaviside step function: 

\ (t>0), 



9{t) = 



(t< 0). 



(18) 



When a wavepacket \ip), containing a pair of photons, passes 
through an atom in the ground state \g), the probability am- 
plitude of the two-photon excitation can be given as 



a 2 = (b,0\U 2 \g,il>) 
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Similar to the case of eq. (1 1), we can define the two-photon 
amplitude as 



ip(h,t 2 ) = (0|o 2 (*2)ai(ti)|^), 



(20) 



whose square corresponds to the joint probability of finding 
photon 1 at t = t\ and photon 2 at t = t 2 . The above two- 
photon amplitude is called an effective two-photon wave- 
function or a biphoton. 17 ' 20 It is clear that the part of 
ip(ti,t 2 ) for t 2 < t\ does not contribute to the two-photon 
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Fig. 2. Schematic diagram of (a) |\P(u>i , and (b) |f (ui , 



excitation because the absorption of photon 1 is always fol- 
lowed by that of photon 2. In addition to the Fourier trans- 
form of tp(ti, t 2 ), which is given as 

1 />oo />OC 
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we also introduce the Fourier transform of ip(ti,t2)0(t2 
ti): 



*(wi,w 2 ) 
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dt 2 dtj 



xe^e^^ti.iaj^-ti); (22) 



this Fourier transform represents the spectral intensity of the 
photon pair under the time ordering i 2 > t\. From eq. (19), 
we derive the probability of two-photon excitation, P 2 = 

|c*2| 2 , as 

P 2 =47r 2 r 1 r 2 |^(0 1 ,f> b )| 2 , (23) 

where r x = £lT>l/h 2 and r 2 = £ 2 V 2 /h 2 . The probability 
of two-photon absorption is given by the spectral component 
|^(f2 a , 51f,)| 2 of the time-ordered wavefunction at the transi- 
tion frequencies of il a and 51^. 

2.4 One-photon transition induced by one photon of the 
photon pair 

In this section, we consider the probability Pi of pho- 
ton 1 exciting the atom to the intermediate level \a). In 
eq. (13), the probability of one-photon absorption is deter- 
mined by the spectrum of the photon that induces the transi- 
tion. If the two photons are not entangled, i. e. , ^(uji, u) 2 ) = 
'J'i(cji) ® ^2(^2), we can easily derive the excitation prob- 
ability as Pi = 27Tr-i 1 (r2 tt ) | 2 . However, in the case of an 
entangled photon, shown in Fig. 2(a), the spectrum of photon 
1 cannot be represented as the function of u>i alone, because 
the spectrum of photon 1 depends on that of photon 2. In 
such cases, the probability of one-photon absorption is ob- 
tained by tracing over the state of photon 2 as 



Pi = 27rn 



(24) 



3. Two-photon absorption induced by time-correlated 
photon pairs 

Before introducing the two-photon wavefuction specific 
to this study, we will discuss the general properties of two- 
photon absorption induced by a time-correlated photon pair. 



The two-photon wavefunction is characterized by two pa- 
rameters: coherent time T and correlation time r. The co- 
herence time T is equal to the coherence time of the pump 
field. The correlation time r (<C T) corresponds to the time 
period required for both the photons to be detected by the 
two detectors and is determined by the spectral width of 
the down-converted light. In the frequency domain, the two- 
photon wavefunction ^(u!%, oj 2 ) can be represented as the 
product of two functions: ^ + (uj + ) with a width of 1/T and 
with a width of 1/r , where u + = u>\ + u>2 and 
w_ = (u>i — lo 2 )/2, as shown in Fig. 2(a). From eq. (22), 
fy(u)i, u> 2 ) is represented as the convolution integral: 

^>(w 1 ,uj 2 ) = / do/_*_(u;_ -c/_)6(o/_), (25) 



where 0(w) is the Fourier integral of the Heaviside step 
function 6(t) and is given by 



e(o;) = </-%)- 1 



2ir w 



(26) 



The symbol P denotes the Cauchy principal value. 21 The 
second term in the above equation causes the spectrum of 
fy(uii, 0J2) to broaden in the direction of the frequency dif- 
ference uj-, as shown in Fig. 2(b), because of its convolution 
integral with 

For simplicity, it is assumed that 

* + (u, + ) = -^=<(>((u, + -cj + )T) 



*_(ct>_) 



ft 



(27) 



where the common function <fr(aoj), which satisfies 
/ (f>(au))du = 1, has a width of 1/a, height of a, and shows 
a peak at u)± . (The coefficient j3 must be determined to sat- 
isfy the normalization condition J duj±\ i B±(uj±)\ 2 = 1.) 
When the two-photon-resonance condition is satisfied as 

lu+ = n+(= n a + n 6 ), we get = PHo)/VT = 

ftVT. If the detuning A is sufficiently large for <j>((oJ- — 
u>-)t) to be regarded as a delta function, we can obtain 
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(28) 



where A = (fl a — flb)/2 — cj_ corresponds to the detun- 
ing from the intermediate state. It should be noted that the 
two-photon absorption probability scales as P 2 oc 1/A 2 . 
This result is consistent with the fact that two-photon pro- 
cesses via a virtual state are a function of 1/A 2 , which is 
derived for two cavity modes. 22 Equation (28) also indicates 
that the time correlation r <C T, which is a unique property 
of time-correlated photon pairs, substantially enhances the 
two-photon absorption probability P 2 . 

Next, we compare the probability of two-photon excita- 
tion, P 2 , with that of one-photon excitation to the intermedi- 
ate state I a), Pi, for the incidence of a photon pair. For any 
two-photon wavepacket, we can calculate P 2 and P x using 
eq. (23) and eq. (24), respectively. The probabilities P x and 
P 2 are significantly affected by the shape of the two-photon 
wavefunction, as described in Sec. 2. In the following sub- 
sections, we investigate these probabilities in two cases: the 
case of a Gaussian wavefunction, as shown in Fig. 3(a), and 
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Fig. 3. Two-photon wavefunction: (a) Gaussian and (b) rectangular. 



the case of a rectangular wavefunction, as shown in Fig. 3(b). 
The Gaussian wavefunction can be obtained by restricting 
the spectrum of photon pairs using a Gaussian filter. The 
rectangular wavefunction can be generated by a photon pairs 
via type-II spontaneous parametric-down conversion due to 
the difference in group velocities of two photons. 15 ' 17 

3.1 Gaussian wavefunction 

The Gaussian two-photon wavefunction is represented as 



x e - iai * 1 e~ i<32t2 , 



(29) 



where Qi (ui 2 ) is the central frequency of the signal (idler) 
photon and a x = 1/\/2tttT is the normalization factor. The 
Fourier transform of this wavefunction is easily derived as 



(30) 



where A\ = ^/2tT/tt. From its definition given in eq. (22), 
we obtain 



w(wi,W2j = — e 



F((u--Q-)t), (31) 



where F(£) is called plasma dispersion function 23 and de- 
fined as 



F(0 



1 



4= /W 



(32) 



It should be noted that the Gauss function in the last factor 
of eq. (30) is replaced with the plasma dispersion function in 
eq.(31). 

By substituting eq. (31) in eq. (23) and assuming that the 
sum of the two-photon frequency is tuned to the two-photon 
transition, i.e., cD + = fl a + fib, we obtain the probability of 
the two-photon transition as 

P 2 = TT 2 r 1 r 2 A 2 1 \F(A-T)\ 2 . (33) 

From eqs. (30) and (23), the probability of the one-photon 
transition to the intermediate state la) can be derived as 
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2(A-r) 2 



T 



(34) 



In the above derivation, we assume that T 3> t, which 
allows us to consider e - 2 (^2-^b) t as ^ delta function 
^j2S{uj 2 -fl b )/T. 
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Fig. 4. (a) Comparison of \F(A ■ r)| 2 (solid line) with c" 2 ( A ' T ) 2 
(dashed line), (b) Comparison of sine 2 (A ■ r/2) (solid line) with 
sine 2 (A ■ r) (dashed line). 



Here, we introduce the ratio 

P 2 fn \F(A-r)\ 2 



(35) 



In the second factor, which is a function of the correlation 
time r, both the Gauss function c~ 2 ' A ' r ) in the numerator 
and the square of the plasma dispersion function |F(A • r) | 2 
in the denominator are monotonously decreasing functions, 
as shown in Fig. 4(a). However, there exists a critical differ- 
ence in the asymptotic behavior of A • r > 1: e~ 2 ( A T ) 
decreases rapidly, while \F(A ■ t)| 2 decreases slowly as 
1/(At) 2 . Hence, we expect that the two-photon transition 
probability P 2 exceeds the one-photon transition probability 
Pi depending on the value of A and selected t, even if the 
first factor in eq. (35) is less than 1. 

3.2 Rectangular wavefunction 

By performing the type-II spontaneous parametric-down 
conversion and using a birefringent crystal to achieve group 
velocity compensation, we prepare the two-photon wave- 
function, which is expressed as 

MhM) = a 2 e-^ +t ^' WT2 

x n T (ti -t 2 )e~ iaitl e - [a2t2 , 



(36) 



where n T (t) is a window function defined as 

'0 (|t|>r), 



n T (t) = 

and the normalization 

lAM^/V^T. 17 - 20 



1 (|t|<T), 

constant is 
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As shown in Sec. 3.1, we can deduce ^(wi,^) and 
^!{uj\,W2) as follows: 

*(wi,w 2 ) = A 2 c- {u1 +- lD + )t2 sinc((u;_ -w_)r) (38) 
*( Wl , W2 ) = ^e-(- + -^) T2 



xsmc| ^ r ) c i(-^W2 



(39) 



where sine a; = sinx/x and A 2 = (2/-k 3 ) 1 / 4 ^/tT. On sub- 
stituting eqs. (38) and (39) in eqs. (23) and (24), respectively, 
we obtain 
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Then, the ratio R r is given as 



R ~ P2 
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The first factor in this equation is identical to that in eq. (35) 
because the two wavefunctions expressed in eqs. (29) and 
(36) exhibit the same function with respect to t\ + t 2 . In 
the second factor, the period of the function of the numer- 
ator is twice as long as that of the denominator, as shown 
in Fig. 4(b). For r = 27m/ A, where n is an integer, the 
two-photon absorption is no longer observed. Fei et al. pre- 
dicted this phenomenon and called it entanglement-induced 
two-photon transparency. 14 From this study, we have found 
that for t = 7r(2n + 1)/A, only the two-photon absorption 
is induced, while the one-photon absorption is completely 
suppressed. 

4. Conclusion 

We derived the probabilities for one-photon and two- 
photon transitions (Pi and P 2 , respectively) for two-photon 
states in general. We showed that the two-photon absorption 
can be dramatically enhanced because of the time correlation 
t « T of general two-photon wavefunctions of a photon 
pair. Then, we dealt with two typical examples of a Gaus- 
sian wavefunction and a rectangular wavefunction and cal- 
culated Pi and P 2 in both cases. The probabilities Pi and 
P 2 were found to behave differently with respect to A • r. 
On the basis of this difference in behavior, we can enhance 



the two-photon absorption while suppressing the undesired 
one-photon absorption, by adjusting the detuning A and the 
correlation time r. In particular, the photon pair having a 
rectangular wavefuction can induce the two-photon transi- 
tion without any one-photon loss under certain special con- 
ditions. 
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